Effects of variable-temperature heat reservoirs on performance of irreversible Carnot refrigerator with heat recovery

The outlet temperature of the heat recovery reservoir is an important parameter in the design of refrigeration with heat recovery systems. In this paper the second law of thermodynamics has been applied to an irreversible Carnot refrigerator with heat recovery (CRHR) coupled to variable-temperature heat reservoirs. The refrigerating rate, input power, refrigeration coefficient, heat recovery coefficient, comprehensive coefficient of performance and exergy efficiency are chosen as the objective functions. The design rule chosen for this study is that the heat transfer area should be constrained. The mathematical expressions for assessing performance parameters with respect to area ratio, were derived for this study. These expressions are transcendental equations. The numerical solution method was employed to calculate the approximate solutions of the optimum performance parameters in a numerical example. The results indicate that the increase in the outlet temperature of heat recovery reservoir could lead to a rise in the maximum value of refrigerating rate and minimum value of input power; also it will lead to the decline in the maximum value of refrigeration coefficient, heat recovery coefficient, comprehensive coefficient and the exergy efficiency. When the ratio of heat recovery heat exchanger area to the summation of high temperature heat exchanger area and the heat recovery heat exchanger area is 1.0, the performance coefficients would attain their limit values and all of the condensing heat could be recycled. Our findings are helpful to the design and optimization to inform preparation of standard relating to the development of refrigerator with heat recovery.


Effects of variable-temperature heat reservoirs on performance of irreversible Carnot refrigerator with heat recovery
Zhe Zhang 1 , Huan Su 1,2* , Guoqiang Dai 1 , Xiaohua Li 1,2 & Liping Zeng 1,2 The outlet temperature of the heat recovery reservoir is an important parameter in the design of refrigeration with heat recovery systems.In this paper the second law of thermodynamics has been applied to an irreversible Carnot refrigerator with heat recovery (CRHR) coupled to variabletemperature heat reservoirs.The refrigerating rate, input power, refrigeration coefficient, heat recovery coefficient, comprehensive coefficient of performance and exergy efficiency are chosen as the objective functions.The design rule chosen for this study is that the heat transfer area should be constrained.The mathematical expressions for assessing performance parameters with respect to area ratio, were derived for this study.These expressions are transcendental equations.The numerical solution method was employed to calculate the approximate solutions of the optimum performance parameters in a numerical example.The results indicate that the increase in the outlet temperature of heat recovery reservoir could lead to a rise in the maximum value of refrigerating rate and minimum value of input power; also it will lead to the decline in the maximum value of refrigeration coefficient, heat recovery coefficient, comprehensive coefficient and the exergy efficiency.When the ratio of heat recovery heat exchanger area to the summation of high temperature heat exchanger area and the heat recovery heat exchanger area is 1.0, the performance coefficients would attain their limit values and all of the condensing heat could be recycled.Our findings are helpful to the design and optimization to inform preparation of standard relating to the development of refrigerator with heat recovery.

List of symbols C
Thermal capacitance rate (kW K −1 ) COP int Comprehensive coefficient F Heat transfer area (m 2 ) f The ratio of the summation of high temperature heat exchanger area and the heat recovery heat exchanger area to the total heat exchanger area f 0 The ratio of heat recovery heat exchanger area to the summation of high temperature heat exchanger area and the heat recovery heat exchanger area K Heat transfer coefficient (kW K −1 m − The exergy efficiency In the last decades, there has been a rapid development to combine space cooling and water heating systems for residences 1 .This combined system is considered to be cost 2 , environmental and energy efficient technology, which recycles some or all of the waste heat in the process of refrigeration for the demand of sanitary water, industrial heating or air reheating and so on 3 .Nowadays, finite-time thermodynamics analysis of thermodynamic systems has become a prominent topic in heat engine, refrigerator and heat pump and so on.Since the 1950s, finite-time thermodynamics has made significant progress after the landmark paper of Novikov 4 and Chambadal 5 .In the 1970s, Finite-time thermodynamics was further advanced by Curzon and Ahlborn 6 .In recent years, many studies on second law of thermodynamics systems have been presented by many authors 7 . Zhang 8 et al. studied the optimization of heat exchanger structure based on finite-time thermodynamics.Li Tao 9 used finite-time thermodynamics to optimize the heat pump system and found that a reasonable selection of the heat transfer area of the heat exchanger can effectively improve the operating performance of the system.John 10 utilized the Carnot cycle to evaluate the potential for waste heat recovery at a wastewater treatment plants.Chen et al. [11][12][13][14][15][16][17] have analyzed the performances of simple and regenerated, endoreversible and irreversible, constant and variable-temperature heat-reservoir air heat-pumps considering the heat load, heat-load density, coefficient of performance (COP) and so on.Sarkar 18 studied the minimization of heat exchanger area or overall conductance of heat pumps and refrigerators for a specified capacity and the analytical results were confirmed by a detailed numerical simulation.Wu 19 proposed an original sinusoidal wavy winglet type vortex generator and evaluated the grade of energy and to explore the irreversible loss during the heat transfer process in view of the second law of thermodynamics.Lei 20 optimized the heat exchanger area by using the finite time thermodynamics theory.
Tan 21 established an endoreversible Carmot cycle model by using finite time thermodynamics.The surface area distribution of three kinds of heat exchangers was optimized by numerical calculation method, and the maximum output power was obtained.Tyagi et al. [22][23][24][25] applied the finite time thermodynamics to various endoreversible and irreversible cycles.They have investigated the effects of a finite rate of heat transfer or other major irreversibilities on the performance of different cycles.Ruibo 26 Applying finite-time thermodynamics theory, an irreversible steady flow Lenoir cycle model with variable-temperature heat reservoirs is established, the expressions of power (P) and efficiency (η) are derived.Based on the theory of finite-time thermodynamics, Meng Fankai 27 designed a channel structure of cooling air, and established a finite-time thermodynamic model of variable temperature heat source thermoelectric cooler based on heat pipe heat dissipation.The thermal resistance of the cold and hot ends of the device was analyzed by numerical simulation method.Wang 28 built an irreversible Carnot heat engine cycle model for space power plants by using finite time thermodynamics.The influences of internal irreversible effect and heat leakage loss on the optimum power output performance are analyzed, when thermal conductivity coefficients of the heat exchanger and cold exchanger are given.Wu 29 studied a reversible simple air refrigeration cycle by using classical thermodynamics.Through theoretical analysis and numerical calculations, the optimal performance of the refrigeration cycle is given.The influence of cycle temperature ratio on the optimal performance of refrigeration cycle is analyzed.The performance of refrigerators with heat recovery coupled to three constant heat reservoirs has been researched using the second law of thermodynamics by our Team 1,5,9,30 .The outlet temperature of the heat recovery reservoir is an important parameter for heat recovery and the variable-temperature heat reservoir is much closer to actual conditions.The performance of the irreversible Carnot refrigerator with heat recovery (CRHR) coupled to variable-temperature heat reservoirs has not been studied by other researchers.El-Din 31 applied the second law of thermodynamics to irreversible heat pumps and refrigerators with two variable temperature heat reservoirs.Q H and Q L were chosen to be the objective functions for heat pumps and refrigerators respectively.However, the results are imperfect due to the optimal variable x (thermal conductance ratio) contained in E h and E c ; these were treated as an invariable during the process of the derivation for maximization of Q H and Q L .In this paper the second law of thermodynamics was applied to an irreversible Carnot refrigerator with heat recovery coupled to variable-temperature heat reservoirs.The refrigerating rate ( R ), input power ( P ), refrigeration coefficient ( ε ), heat recovery coefficient ( ε R ), comprehensive coefficient of performance ( COP int ) and exergy efficiency ( η � ) were chosen as the objective functions in this study.Equating the derivatives of those performance parameters with respect to area ratio, f, to zero, a group of transcendental equations would be derived.The numerical solution method was employed to calculate the approximate solutions of the optimum performance parameters in a numerical example.The influence of outlet temperature of the heat recovery reservoir on these performance parameters was analyzed in the numerical example.

Thermodynamic model
An irreversible Carnot refrigerator with heat recovery coupled to variable-temperature heat reservoirs and its surroundings are shown in Fig. 1a.There are three heat exchangers, including high temperature heat exchanger, low temperature heat exchanger and heat recovery heat exchanger, existing in the cycle.
Figure 1b shows the temperature-entropy diagram of an irreversible Carnot refrigerator with heat recovery which operates steadily between three variable temperature reservoirs.The working fluid in the refrigerator has two constant temperatures, T HC and T LC .The three heat exchangers are treated as counter flow heat exchanger with finite thermal capacitance rates,C r ,C h and C l .
In Fig. 1, A-the compressor; B-the heat recovery heat exchanger; C-the high temperature heat exchanger; D-the expander, E-the low temperature heat exchanger.
Generally, for a combined space cooling and water heating system, the inlet temperatures of high temperature reservoir ( T Hi ) and heat recovery reservoir ( T Ri ) can be treated as the outdoor environment temperature as described in Eq. ( 1) and the inlet temperature of low temperature reservoir ( T Li ) can be treated as the indoor environment temperature as given in Eq. ( 2).
(1) where T o and T i are the outdoor and indoor environment temperature.
The rate of heat transfer at the high temperature side (Q H ) is defined by Eq. ( 3): From Eq. (1), Eq. (3) becomes Eq. ( 4): where C h is the thermal capacitance rate of high temperature reservoir and T Ho is the outlet temperature of the high temperature reservoir.
Using the LMTD method, Q H can also be determined by Eq. ( 5): where U h is the heat transfer coefficient of high temperature reservoir.F h is the heat transfer area of high tem- perature reservoir.T HC is the temperature of working fluid at the high temperature side.From Eqs. ( 4) and ( 5), Eq. ( 6) is derived.
From Eqs. ( 6), ( 7) is derived for determination of the outlet temperature of the high temperature reservoir.
where η h is the effectiveness of the high temperature heat exchanger; and is defined by Eq. ( 8): Substituting Eq. ( 7) into Eq.( 4) gives Eq. ( 9): The rate of heat transfer at the heat recovery side (Q R ) is determined by Eq. ( 10): Equation ( 11) is derived from Eq. (1) and Eq. ( 10) for determination of Q R : where C r is the thermal capacitance rate of heat recovery reservoir.Using the LMTD method, Q R can also be determined by Eq. ( 12), as: where U r is the heat transfer coefficient of heat recovery reservoir.F r is the heat transfer area of heat recovery reservoir.From Eqs. (11) and 12), Eq. ( 13) is derived: From Eqs. ( 13), ( 14) is derived: where η r is the effectiveness of the heat recovery heat exchanger; and is defined by Eq. ( 15).
Substituting Eq. ( 14) into Eq.(11) gives Eq. ( 16): The rate of heat transfer at the low temperature side, Q L is determined by Eq. ( 17): From Eqs. ( 2), ( 18) can be obtained to determine Q L : where C l is the thermal capacitance rate of low temperature reservoir.Using the LMTD method, Q L can also be determined by Eq. ( 19), as: (2) where U l is the heat transfer coefficient of low temperature reservoir.F l is the heat transfer area of low temperature reservoir.From Eqs. ( 18) and ( 19), Eq. ( 20) is derived: From Eqs. ( 20), ( 21) can be derived: where η l is the effectiveness of the low temperature heat exchanger; and is defined by Eq. ( 22).
Substituting Eq. ( 21) into Eq.( 18) gives Eq. ( 23): The cycle is internally irreversible and the internal irreversibility parameter can be determined by Eq. ( 24): where is the internal irreversibility parameter, which is always greater than 1 for an irreversibility cycle and equals to 1 for an endoreversible cycle.
The ratio of the heat transfer rate at the heat recovery side to the total heat emissions, n can be defined by Eq. ( 25), below: where, n is defined as the heat recovery ratio.
The design rule chosen by this paper is that the heat transfer area should be constrained, as defined by Eq. ( 26): The ratio of F h and F r to F is defined by Eq. ( 27): And Eq. ( 28) defines the ratio of F r to F h and F r .
where f is the ratio of the summation of high temperature heat exchanger area and the heat recovery heat exchanger area to the total heat exchanger area.f 0 is the ratio of heat recovery heat exchanger area to the sum- mation of high temperature heat exchanger area and the heat recovery heat exchanger area.From Eqs. ( 26), (27) and Eqs. ( 28), ( 29), (30) and (31) are respectively derived: The equations above constitute the mathematical model for optimizing the performance of an irreversible Carnot refrigerator with heat recovery operating between the three variable-temperature heat reservoirs.From these equations, Eqs.(32) to (36) can be derived for description of the objective functions as given below: www.nature.com/scientificreports/I n t h e s e e q u a t i o n s , Where R is the refrigeration rate; P is the input power; ε is the refrigeration coefficient; ε R is the heat recovery coefficient; which is defined as the ratio of heat recovery rate to the input power.COP int is the comprehensive coefficient, which is the summation of refrigeration coefficient and heat recovery coefficient.E p is the input elec- tricity exergy of the CRHR; E L is the cold exergy; E R is the recovery heat exergy; E H is the heat exergy emitted to the surroundings directly, which is a part of the exergy loss; is the other exergy loss of the CRHR.
In Eqs.(32) to (37) the superscript point on f means that f is chosen to be the optimization variable when the parameters such as C h ,C r ,C l ,U h ,U r ,U l ,F,T i ,T o ,T Ro , f o and are specified.Maximizations or minimizations of these performance parameters with respect to f give Eqs.(38) to (43): By calculating dR/df = 0,dP/df = 0,dε/df = 0,dε R /df = 0 ,dCOP int /df = 0 and dη � /df ≥ 0 ; the maxi- mums/minimums of R,P,ε,ε R , COP int and η � can be obtained.From Eqs. (32) to (43), six of the objective func- tions are obviously the transcendental equations and the analytical formulas cannot be obtained for the optimum performance parameters.In this paper, the numerical solution method was used to calculate the approximate solution by Matlab software.The influence of temperature variations of heat recovery reservoir on the optimal performance parameters can be researched by substituting different values of T Ro .In this way, the influence of f 0 and also can be studied.

Results and discussion
The specified parameters and numerical examples are listed in Table 1 30,31 : There are two main input variables in this paper, one of which is that the value range of the ratio of the summation of high temperature heat exchanger area and the heat recovery heat exchanger area to the total heat exchanger area ( f ) is 0-1, and the other is T Ro , which is the outlet temperature of heat recovery reservoir.Because the heat recovery in this paper is used to produce sanitary water and industrial heating.Therefore, the lowest temperature of the T Ro must be higher than the ambient temperature, and the highest temperature should not exceed the temperature of the condenser in the refrigeration cycle.Hence, in this paper, the temperature range of the T Ro is selected as 310-370 K. www.nature.com/scientificreports/

Optimal value of R
The plot of R vs. f is shown in Fig. 2 for three values of T Ro under the condition that: = 1 and f 0 = 0.5.The plot of f −opt−R vs.T Ro is shown in Fig. 3 for three values of f 0 and two values of .In Fig. 2, the optimal values of f denoted as f −opt−R , at which the refrigerating rate would attain their maxi- mum values, are denoted as R f .In this case the optimal values of f are f −opt−R = 0.508, 0.505 and 0.503; the maximums of refrigeration rate are R f = 80.07 kW, 103.45 kW and 121.18 kW corresponding to T Ro = 328 K, 338 K and 348 K respectively.The greater is the T Ro , the greater would be the refrigeration rate.Figure 3 shows that T Ro , and f 0 all have effects on f −opt−R .The value of f −opt−R moves closer to a constant value gradually, along with an increase in T Ro .This constant value is influenced by .The smaller is the value of f 0 , the smaller would be the f −opt−R .
The plot of R f vs.T Ro is shown in Fig. 4 for three values of f 0 and two values of .Figure 4 shows an increase in R f with an almost corresponding increase in T Ro .The smaller is the value of f 0 , the greater would be the R f .The greater is the value of , the smaller would be R f .

Optimal value of P
The plot of P vs. f is shown in Fig. 5 for three values of T Ro under the condition that: = 1 and f 0 = 0.5.In Fig. 5, the optimal values of f denoted as f −opt−P , at which the input power attain minimum values are denoted as P f .In this case the optimal values of f are f −opt−P = 0.467, 0.465 and 0.463.The minimum values of input power are P f = 43.86kW, 82.38 kW and 126.58 kW, corresponding to T Ro = 328 K, 338 K and 348 K respectively.The greater is the T Ro , the greater would be the input power.
The plot of f −opt−P vs. T Ro is shown in Fig. 6 for three values of f 0 and two values of .Figure 6 shows that T Ro , and f 0 all have effects on f −opt−P .Different from the curves of f −opt−R vs.T Ro , the curves of f −opt−P diverge from each other along with an increase in T Ro .www.nature.com/scientificreports/ The plot of P f vs. T Ro is shown in Fig. 7 for three values of f 0 and two values of .Figure 7 shows that P f increase along with an increase in T Ro .The smaller is the value of f 0 , the greater would be the P f .The greater is the value of , the greater is P f .

Optimal value of ε
The plot of ε vs. f is shown in Fig. 8 for three values of T Ro under the condition that: = 1 and f 0 = 0.5.In Fig. 8, there is only one optimal value of f , denoted as f −opt−ε , at which the refrigeration coefficient attains maximum values, denoted as ε f .In this case the optimal value of f is f −opt−ε = 0.483; the maximum values of refrigeration coefficient are ε f = 1.83, 1.26 and 0.96 corresponding to T Ro = 328 K, 338 K and 348 K respectively.The greater is the T Ro , the smaller would be the refrigeration coefficient.
The influence factors on f −opt−ε can be analyzed by the following plots.Figure 9a shows the plot of f −opt−ε vs. T Ro for three values of f 0 and two values of .T Ro is shown to have no effect on f −opt−ε .Figure 9b shows the plot of f −opt−ε vs. f 0 under the condition that: = 1 and T Ro = 328 K.In the range of the available value of f 0 ,0 ≤ f 0 ≤ 1 , f −opt−ε has a maximum value and a minimum value.When f 0 = 0, f −opt−ε attains its maximum value (about 0.55); when f 0 = 0.7, f −opt−ε attains its minimum value (about 0.481).Figure 9c shows the plot of f −opt−ε vs. under the condition that f 0 = 0.5 and T Ro = 328 K. f −opt−ε vs. is approximate to the monotonous linear relationship.The greater is the , the greater would be the value of f −opt−ε .
The influence factors on ε f can be analyzed by the following plots.Figure 10a shows the plot of ε f vs. T Ro for three values of f 0 and two values of .
The maximum ε f decreases dramatically along with the increase of T Ro .From Eq. ( 11) the heat recovery rate Q R is shown to increase monotonically with T Ro due to the specified values of T o and C r .That means the greater is T Ro , the larger would be the recycled heat and the lower is refrigeration coefficient.Figure 10b shows the plot of ε f vs. f 0 under the condition that = 1 and T Ro = 328 K and Fig. 10c shows the plot of ε f vs. under the condition www.nature.com/scientificreports/ that f 0 = 0.5 and T Ro = 328 K. ε f vs. f 0 as well as ε f vs. are approximate to the monotonous linear relationships.The greater is the f 0 , the greater would be the value of ε f .The greater is the , the smaller would be ε f .

Optimal value of ε R
The plot of ε R vs. f is shown in Fig. 11 for three values of T Ro under the condition that: = 1 and f 0 = 0.5.In Fig. 11, the optimal value of f is denoted as f −opt−εR , at which the heat recovery coefficient attains its maximum values, denoted as ε Rf .In this case the optimal values of f are f −opt−εR = 0.467, 0.465 and 0.463.The maximum values of refrigeration coefficient are ε f = 1.37, 1.09 and 0.95 corresponding to T Ro = 328 K, 338 K and 348 K respectively.The values of f −opt−εR is the same as f −opt−P due to the definition of ε R in Eq. ( 35).The greater is the T Ro , the smaller would be the heat recovery coefficient.
The influence factors on f −opt−εR can be analyzed by the following plots.Figure 12a shows the plot of f −opt−εR vs.T Ro for three values of f 0 and two values of .The curves in Fig. 12a are the same as Fig. 6. Figure 12b shows the plot of f −opt−εR vs. f 0 under the condition that = 1 and T Ro = 328 K.In the range of the available value of f 0 , ( 0 ≤ f 0 ≤ 1) f εR has a minimum value.When f 0 = 0, f −opt−εR gradually tend to be 1; when f 0 = 0.72, f −opt−εR attains its minimum value (about 0.46).Figure 12c shows the plot of f ε R vs. under the condition that f 0 = 0.5 and T Ro = 328 K. f −opt−ε vs. is approximate to the monotonous linear relationship.The greater is the , the greater would be the value of f −opt−εR .
The influence factors on ε Rf can be analyzed by the following plots.Figure 13a shows the plot of ε Rf vs. T Ro for three values of f 0 and two values of .The maximum ε Rf declines sharply along with an increase in T Ro , due to the input power P is increasing faster than Q R .Figure 13b shows the plot of ε Rf vs. f 0 under the condition that = 1 and T Ro = 328 K and Fig. 13c shows the plot of ε Rf vs. under the condition that f 0 = 0.5 and T Ro = 328 K. ε f vs. f 0 as well as ε f vs. are approximate to the monotonous linear relationships.The influence rules of f 0 and on ε Rf is similar to ε f .

Optimal value of COP int
From definition of COP int in Eq. ( 36) that there must be an optimal value of f , denoted as f −opt−COPint , at which the comprehensive coefficient could attain its maximum values, denoted as COP intf .attains its minimum value (about 0.473).Figure 14c shows the plot of f −opt−COPint vs. under the condition that f 0 = 0.5 and T Ro = 328 K. f −opt−COPint vs. is approximate to the monotonous linear relationship.The greater is the , the greater would be the value of f −opt−COPint .
The influence factors on COP intf can be analyzed by the following plots.Figure 15a shows the plot of COP intf vs. T Ro for three values of f 0 and two values of .The maximum COP intf declines sharply along with an increase in T Ro .Figure 15b shows the plot of COP intf vs. f 0 under the condition that = 1 and T Ro = 328 K and Fig. 15c shows the plot of COP intf vs. under the condition that f 0 = 0.5 and T Ro = 328 K. COP intf vs. f 0 as well as COP intf vs. are approximate to the monotonous linear relationships.Figures 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, and 15 show the completed approximate solutions of a numerical example.The results of this numerical solution example are under the specified conditions, but the influence factors on the performance parameters are of general significance.With specified heat transfer coefficients and thermal capacitance rates there must be a maximum value of the refrigerating rate, refrigeration coefficient, heat recovery coefficient or comprehensive coefficient and a minimum value of the input power existing in the cycle.The increase of the outlet temperature of heat recovery reservoir could lead to a rise in the maximum value of refrigeration rate and minimum value of input power as well as to the decline of the maximum value of refrigeration coefficient, heat recovery coefficient and comprehensive coefficient.The rise of f 0 is beneficial to the performance coefficients.The rise of is harmful to them.

Optimal value of η �
The plot of η � vs. f is shown in Fig. 16 for three values of T Ro under the condition that: = 1 and f 0 = 0.5.The plot of f −opt−η � vs.T Ro is shown in Fig. 17 for three values of f 0 and two values of .
In Fig. 16, the optimal values of η � denoted as f −opt−η � , at which the exergy efficiency would attain their maximum values, are denoted as η � f .In this case the optimal values of η � are f −opt−η � = 0.4942, 0.4958 and 0.4927; the maximums of exergy efficiency are η � f = 0.2035, 0.1532 and 0.1248 corresponding to T Ro = 328 K, 338 K and 348 K respectively.The greater is the T Ro , the smaller would be the exergy efficiency.Figure 17 shows that T Ro , and f 0 all have effects on f −opt−η � .The value of f −opt−η � moves closer to a constant value gradually, along with an increase in T Ro .This constant value is influenced by .The smaller is the value of f 0 , the smaller would be the f −opt−η � .

Conclusion
In this paper, we applied the second law of thermodynamics to an irreversible Carnot refrigerator with heat recovery coupled to variable-temperature heat reservoirs.The heat recovery process involves the recycling of waste heat generated during refrigeration for various purposes, such as sanitary water supply, industrial heating, and air condensation.Through the use of numerical solutions, we were able to obtain accurate results for a specific example.Our study yielded several key findings, which are summarized below: (1) With specified heat transfer coefficients and thermal capacitance rates there must be an optimal value of f at which the performance parameter attains maximum or minimum value.The outlet temperature of heat recovery reservoir ( T Ro ) can have an effect on f −opt−R , f −opt−P , f −opt−ε R , f −opt−COP int and f −opt−η � ; but it has no effect on f −opt−ε .(2) The increase in the outlet temperature of heat recovery reservoir ( T Ro ) could lead to a rise in the maximum value of refrigerating rate ( R f ) and minimum value of input power ( P f ); also it will lead to the decline in the maximum value of refrigeration coefficient ( ε f ), heat recovery coefficient ( ε Rf ), comprehensive coefficient ( COP intf ) and the exergy efficiency ( η �f ).(3) The rise of f 0 is beneficial to the performance coefficients, but it could lead to a decline in R f .When f 0 = 1.0 the performance coefficients would attain their limit values and all of the condensing heat could be recycled.The rise of can be harmful to the performance coefficients.

Figure 1 .
Figure 1.(a) The schematic of the refrigerator with heat recovery, (b) The temperature-entropy diagram of the CRHR cycle.

Figure 2 .
Figure 2. R vs. f for three values of T Ro .

Figure 3 .Figure 4 .Figure 5 .
Figure 3. f −opt−R vs. T Ro for three values of f 0 and two values of .

Figure 6 .
Figure 6.f −opt−P vs. T Ro for three values of f 0 and two values of .

Figure 7 .Figure 8 .
Figure 7. P f vs. T Ro for three values of f 0 and two values of .

Figure 16 .
Figure 16.η � vs. f for three values of T Ro .